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SUMMARY 

Using the lifting surface theory and the acceleration potential method for 
the flow field of an axial turbo-compressor stage, a recursive and a direct 
method are presented that make use of the eigenfunction solutions of the Iso- 
lated rotor and stator to solve for the rotor-stator Interaction problem. The 
net pressure distribution on the rotor and stator blades Is represented by mod- 
ified Blrnbaum series, whose coefficients are determined using a matrix proce- 
dure and satisfying the boundary conditions on the surface of the blades. The 
relation between the matrix operators of the recursive and the direct methods 
Is also shown. Expressions have been given for the blade circulation, the 
axial and tangential forces on the blade, the rotor power required, and the 
Induced upwash velocity of the stage. 


1 . INTRODUCTION 

The performance estimation of a turbomachine stage requires the aerody- 
namic characteristics of the rotor and the stator including their mutual 
Interaction. The purpose of this paper Is to discuss the aerodynamics of a 
rotor-stator combination In a turbomachine. The spacing between the adjacent 
blade rows of a turbomachine Is generally small compared to the blade chord. 
Hence, the aerodynamic Interference between the rotor and stator blade rows may 
be expected to be significant and needs to be Included for an accurate analysis 
of the stage aerodynamic behavior. Existing theories like those of McCune 
(ref. 1), Namba (ref. 2), Schulten (ref. 3) and others consider only a single 
annular row of blades while other theories like that of Kajl and Okazaki 
(ref. 4) consider two adjacent two-dimensional Infinite cascades. Hence, a 
method Is required for calculating the aerodynamic flow field due to the the 
simultaneous presence of two rows of closely spaced annular cascades and Is 
outlined In the following. 

Since the chordwlse and radial pressure distribution on the rotor and 
stator blades can be represented by a Blrnbaum series, the heart of the lifting 


This paper Is declared a work of ihe U.S. Government and Is 
nol subject to copyright protection In Ihe United States. 



surface method Is the determination of the coefficients of these series. In 
reference 5, these coefficients were determined directly by satisfying flow 
tangency conditions on the upper and lower surfaces of the rotor and stator 
blades simultaneously. It was possible to obtain the parameters governing the 
blade loading and the discrete frequency noise field at the operating point. 

In the following, a recursive scheme for determining the flow field by 
combining the solutions of the Isolated rotor and stator will be outlined. 
Consequently, the procedure for calculating the off-design performance of a 
given rotor-stator combination Is greatly simplified. The recursive procedure 
described here enables one to determine the Blrnbaum coefficients of the rotor- 
stator combination to any desired order of approximation. It Is possible to 
estimate the aerodynamic Interference between the two blade rows by this 
procedure. 

A simple and direct method which solves for the flow field of the rotor 
and stator simultaneously Is also Indicated. 

In section 2, a brief outline of the results of the acceleration potential 
method for dealing with the rotor-stator problem Is given. In section 3, the 
basic solutions to the rotor and stator are obtained by solving a nonhomo- 
geneous problem giving the pressure distribution over the respective blades 
Ignoring aerodynamic Interference effects. These basic coefficients are modi- 
fied to correct for the Interference effects Iteratively. The recurrence pro- 
cedure used to calculate the Blrnbaum series coefficients by satisfying the 
surface boundary conditions on the blades Is discussed In section 3 which also 
contains a direct method for the same purpose. The estimation of the rotor- 
stator aerodynamic Interference Is contained In section 4. In section 5, a dis- 
cussion of the equations for obtaining the circulation distribution over the 
blades, the axial and tangential forces and their distribution, the shaft power 
absorbed by the rotor and the upwash field of the stage. 


2. EIGENFUNCTIONS OF THE ROTOR AND STATOR 

In this section, a brief procedure for obtaining the flow field of the 
rotor and stator using a distribution of acoustic monopole and dipole flow sin- 
gularities on the blade surface Is outlined. The rotor and stator with 
and blades, respectively, are considered to be situated In an unsteady, 

Invlscld compressible flow field with the rotor placed upstream of the stator. 
The Incoming flows for the rotor and stator are assumed to vary harmonically 
with the respective frequencies w r and w s . The perturbations produced by 
the stage are considered to be small compared with the mean flow W a so that 
the linearized equations of motion In three dimensions may be used to represent 
the rotor and stator blades using suitable acoustic singularities. A coordinate 
system rotating with the rotor blades Is chosen with respect to which the free 
stream has the mean velocity components 

(0, V r , W r ) = (0, ar, W a ) (2.1) 

In the cylindrical coordinate system. The axial velocity W a Is assumed to 
be constant through the stage. Relative to the rotating coordinate system 
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fixed to the rotor, the Inlet condition to the stator corresponds approxi- 
mately to the value obtained from the stage velocity diagram and has the 
components 


(0, V s , W s ) > (0, Or - tan a 2 r , W a ) (2.2) 

Since the chordwlse and radial pressure distribution on the rotor and the 
stator blades Is represented by a Blrnbaum series, the heart of the lifting 
surface method Is the determination of the coefficients of this series. In an 
earlier paper, the linearized aerodynamics of the combined rotor-stator system 
of an axial turbomachine for a uniform free stream using an acceleration poten- 
tial representation was determined directly. A Green's function was Introduced 
to satisfy the linearized differential equation for the rotor and stator. The 
radial variation of the Green's function for the rotor follows a Bessel's dif- 
ferential equation. For the stator, due to the presence of a constant swirl 
term W a tan « 2 r Introduced by the rotor due to Its exit blade angle a 2 r , the 
radial variation of the Green's function Is governed by a Whittaker differential 
equation. To make use of convenient orthogonality properties, this equation Is 
transformed Into the Laguerre differential equation. The Bessel and the gener- 
alized Laguerre and Kummer functions are used as the orthonormal basis func- 
tions to form the rotor and stator pressure eigenfunctions for the subsequent 
lifting surface theory. The rotor and stator blades are represented by a sur- 
face distribution of acoustic dipoles. The normalized radial pressure eigen- 
function of the pulsating acoustic unit pressure monopoles on the rotor and 
stator may be shown (ref. 5) to be 

•kfcOV - A(k,ft) J k (x kJl r 1 ) ♦ B( k , a) Y k (x |cJl r 1 ) 

*lkl (r 1 ) = C(l,k,l)L£ 2a) (r 2 ) ♦ D r^F^-l - 25; 1 - 2 i; r 2 ), r 2 = 


(2.3) 

where $kl( r l) the normalized cylinder function In which J k (x), Y^( x) are 
Bessel functions of the first and second kind, respectively; ’Pfckjt(ri) Is the 
normalized stator eigenfunction In which l£“J(x) and -j F-| ( a , b ; x) are, respec- 
tively, the generalized Laguerre function and the Kummer function; and 
v Ikk are the corresponding eigenvalues to be determined; A, B, C, 0 are con- 
stants; k = 0,±1 ,±2 ; 4 = 1,2,3,...; 4 = 0,1, 2, 3 The eigenvalues 

and v lk4 are obtained by satisfying the Neumann boundary conditions 

d§ k o/dr-| =0 at r-\ = h r and 1 

(2.4) 

dY^ki/d i"i ~ ^ 4t f"i = R sr^s 3nd R sr 

at both the hub and tip of the rotor and stator. The pressure field of the 
surface distribution of pulsating unit pressure dipoles on the Isolated rotor 
and stator blades Is obtained by differentiating the monopole pressure function 
normal to the surface. The resulting pressure field for both the rotor and 
stator can be shown to be given by 
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(2.5) 


P r<VV - 



AP r §*(r 1 ,P r ,t 1 ;k r ,i)dp r dc r 


P s (r rV 


f-'f” 

8irMM ~k H I AP s ’M r rV t l ; * > 


k s ,ii)dp s dc s 


where AP r and AP S are the net pressure distribution on the rotor and 
stator blades, respectively which are functions of (p r , <pr» Cr) and 
(ps* «ps. Cs)* The functions $* and V* are defined by 


♦*(r 1 ,P r ,t 1 ;k r ,il) = # k a (p r )^ 0 (r l )f 2 exp (1f o 5 r ) exp 1[k r (e " *r “ v r )J 


( 2 . 6 ) 


**(?!, p s . tjjl.k^i) = ^ p s^0^ r l )f 5 exp <~W> 6XP 1[k s (e " ‘ 


for subsonic axial flow, M < 1. The pressure field P r of the rotor has been 
obtained by summation over all the % x rotor blades and likewise, the field of 
the stator has been obtained by summation over all the stator^blades . The 

functions fo, f2. -^20 and f 5. *^50 are listed In the Appendix; 2T r = z-| - £ r , 
^s = Z 1 " Cs* Pr(pr*<P r .Cr) and Ps(Ps»*s»C$) are the respective locations of 
the dipoles on the rotor and stator blades. 

The radial and chordwlse distribution of pressure on a blade of the rotor 
and the stator Is assumed to be given by the Glauert-BIrnbaum series modulated 
radially by the corresponding radial eigenfunction. The Blrnbaum series coef- 
ficients are assumed to be Independent of the eigen-numbers X, k, it. Assuming 
that the thickness and camber effects are linearly superposable within the 
limits of the linear theory used here, the net pressure distribution on the 
rotor and stator blades can be written as 


AP 


AP 


r 


s 


H 1 (w,P r ;k r ,t) = (ja£ cot | * 
H 2 («,P s ;l,k s ,t) = (& Q cot | 


00 

J sf sin M« k ,(o r > 
m= 1 r 

00 

. I 0 Sin nC)* ikl („ s ) 
m=l 


(2.7) 


4 



These are written In the matrix form 


where 


4 (p r ) 

H 2 («,P s ;i,k s ,i) = t (p s ) 

n, sf,3i denote the vectors 


( 2 . 8 ) 


n = |cot | , sin Z, sin 2«, .... sin mZ, . ..j 

j* T = Cafc.J*,. (2.9) 

# T = (^j,^, & 2 , ...) 

The superscript T Indicates the transpose of the matrix and w Is the 
Glauert angle parameter for the chordwlse position y' ( In the blade coordinate 
system shown In (fig. 1). The coordinates yj. and y^ of the rotor and 
stator are defined by 


y r * -C R cos 


-C 


< y < 

r - J r - 


< ♦ C, 


0 < « r < * 


y s . tc s cos » s 


-c. < y < * c, 


0 < Uj < f 


( 2 . 10 ) 


In order to satisfy the boundary conditions on the blade surfaces, It Is 
necessary to calculate the resultant velocity at a point on the blade surface. 
This Is done by using the unsteady equation of motion In a helical coordinate 
system (fig. 2) and Integrating the perturbation pressure along the undisturbed 
streamline direction (ref. 2) assuming that the pressure and the associated 
velocities also vary harmonically In time. 


3. THE BASIC SOLUTION 

The Blrnbaum coefficient vectors Introduced In section 2 are unknown and 
have to be determined by satisfying the appropriate conditions on the blade 
surfaces. The procedure for calculating the resultant velocity at any point in 
the flow field of the rotor and stator Is described In this section. From this 
the blade surface conditions can be obtained. We shall denote s/(Q) and ^(0) 
to be the basic Blrnbaum coefficient vectors of the Isolated rotor and stator, 
so that mutual Interference effects are absent. In subsonic flow, these coef- 
ficients will be perturbed by aerodynamic Interference effects when the rotor 
and stator are juxtaposed to form a stage. The resulting Blrnbaum coefficients 
of rotor-stator combination will be denoted by si and . The basic coeffi- 
cients j^(°) and are determined by satisfying the surface boundary con- 

ditions on the rotor and stator blades. As the stator approaches the rotor, the 
Induced velocity due to the stator disturbs the surface boundary conditions on 
the rotor. Hence, the rotor Blrnbaum coefficients s/(°) must be changed to 
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satisfy the surface conditions again. Likewise, the stator Blrnbaum coeffi- 
cients £j?(0) must be changed due to the rotor Interference field. 


When the rotor and stator blades are kept In Isolation, the resultant 
velocity 3 r and 3s the blade surface can written In the respective blade 
coordinate systems as 


^rr 


/ u «* A 


sin e 

cos e 


0 

U Ry' 

= 

sin a cos e - 
r 

sin a sin 
r 

6 

COS a 



-cos o p cos e 

cos <* r sin 

6 

sin a p 

X 

(/) 


sin e 

cos e 


0 “ 

U Sy ' 

= 

sin a cos e - 
r 

sin a sin 
r 

e 

COS a 

S 

Wi 


-COS a $ COS 0 

cos o $ sin 

e 

sin 

x/ 

*fe ' (x 2 cos 6 hl 

- X 3 s1n ®M 



(Hr 



*Ve 

+ cot 

6hl , 


+ 1 


/Hr 


\ 

<D 

+ cot 

«h2 


\Hz + A rs J 


( 3 . 1 ) 


Hr 



U, = (X, sin e. 


*se * (Y 2 tos 6 h 2 - Y 3 s1n V* 8 


* X, 


cos e 



Hz ■ (y 2 s1 " ®h 2 * V 3 cos 

where s 4 and OR are the vectors of Blrnbaum coefficients and the functions 
X-| , X2, X3, Yi , Y2, Y3 are listed below. 
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“k»<'r>W<r 


(, - K r I I i' Z ( 

HA 

<2 ' "r l j[| "'I I »k,<'r>-W<r 


' i si "l 




y, . K s 


III f 52 f 

t k 1 1 1 . 


rl 

R 

sr 


(3.2) 


m L^2 e ^ ^ s dC s 


Z sl ' R sr h s 


K .III 

l k l 


f Z s2 r R sr 
%1 *Vs 


2 i “V 2 
"W's^ *5 d 's d *s 


Y 3 = K s 


n sr 

h 

si sr n s 


m, !k»^s)' , 2 e ’ P2> Wp s <tf s 


The functions 9 3, ^4, ^5, ^5, K r , K s are defined In the Appen- 

dix. Figure A1 shows the angles ^r* <Pr for the rotor and likewise Tp s , <p s 
on the stator. If Ri(n, e^, z^) and S^(r^, e^, z^) denote the vector posi- 
tions of a set of points on the rotor and stator, the respective Blrnbaum coef- 
ficients are determined by satisfying the boundary conditions at these points 
on both the upper and lower surfaces. If z‘ = z'(y') Is the equation of the 
blade profile In the blade coordinate system, the flow tangency condition Is 
expressed as t * (dz'/dy 1 ) = ( U z * /Uy 1 ) . Denoting by = 0+ and = 0-, 
the upper and lower surface chord lines, the boundary conditions for the rotor 
and the stator may be written as 


r U 


Rz' 


'U 


Sz' 


- T 


- T 


U1 

U Ry' 

)z'= 0 + 

= 0 

V) 2 ..o- * 0 

U 2 

U Sy' 

)z'= 0 + 

* 0 

(“**• ■ V) 2 , =0 . ■ 0 


(3.3) 


where , x L1 are the chordwlse slopes of the upper and lower surfaces at a 

point on a rotor blade; x u2 and t L2 represent the corresponding quantities 

at points on a stator blade. Substituting from equation (3.1) we can write 
these equations as 
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(3.4) 


Lj/ (0) =^2 U^ (0) =^3 V^<°> = 

which can be rewritten as 



The matrices K, L, U, V, , T 2 are defined In the Appendix. From equa- 
tion (3.5) we obtain the basic Blrnbaum coefficients of the Isolated rotor and 
stator as: 

^ (0) = K” 1 r 1 68 (0) = U v V 2 (3.6) 

It Is possible to calculate the aerodynamic parameters of the rotor and 
stator flow field using the basic coefficients. 


4. THE RECURSION PROCEDURE 

The Blrnbaum coefficients of the Isolated rotor and stator were determined 
In section 3. It Is possible to calculate the Blrnbaum coefficients of the 
rotor- stator combination using a recurrence procedure to be described In this 
section. When the rotor and stator are juxtaposed, the perturbation velocities 
at the rotor and stator blades may be written 
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A w\ if 


Ry' 


V 


= ©_ 


Rz ' / 


,(R) 


\4 R) sin si * *r «« *„ ,/ 
/ *< s > \ 


,<R) 


cos e u , - X< R > sin e, 


\ 


hi 


hi 


(R) 


( 


,(») 


IA T . 


4 R) cos ®h2 - V 3 R) s1n ®h2 ^ * 


/». A 


Sx 


sy 


V 


Sz 


•/ 


= ©. 


( (S) 

'2 

,<S) 


cos 6 ., - xi S) sin 0 , 


hi 


hi 




V 2 


(S) 


sin 0 hl + Xg cos 0 


( 


2 R> s1 " °h2 * ¥ 3 R) cos 6 h2/ 


,(S) 


-y 


\ 


v< 5) cos 0„ 2 - v< s) sin e h2 


,Y< S > sin e h2 *■ cos e ( 


\ 


e T ♦ 


« / 


cot © hl 

V 1 ) 

M 

cot 0^ 2 

W 7 


e r " 


sin 0 


sin a cos © 
r 


,-cos <* r cos 0 


e 


sin 0 


sin <* s cos 0 


^-cos a $ cos 0 


cos 0 


-sin a sin 0 
r 


cos <* r sin 0 


cos 0 


-sin <* s sin 0 


cos a $ sin 0 



(4.1) 


where and Q& are the modified Blrnbaum coefficients written In series form 
as 


= j/(°) + * j^( 2 ) ♦ ... 

(4.2) 

0 = gg(0) + 0(1) * 0(2) * ... 

Using the equations (4.1) and (4.2) In the boundary conditions 
equation (3.4), we obtain 

K(j/(0) ♦j/O) +j/ 2 ) ♦ ...) * S(«(0) t^d) + #(2) * ...) = 

L(^(0) ♦ J /(l) +J /2) * ...) ♦ T(#(0) *-#0) *$0f(2) + ...) =<<f 2 

(4.3) 

H(j*(0) *^(1) +jf( 2) * ...) ♦ U(#(0) ♦#(!) *#<2) + ...) =«gf 3 

N( jtf(O) t^d) +jj^( 2 ) ♦ ...) ♦ V( 0 (O) + #(D + #< 2 ) ♦ ...) =^ 4 

ft 

The four equations In equation (4.3) may be written conveniently as a pair 
of equations 
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K L (a/ 0) +^ (1) +J 2) * ♦ S T ^f 0) +&f 2) + ...) = r 1 


+j* 0) *J 2) * ...) * u ¥ tf> +ef }) *# (2) * ...) = r. 


From equation (3.5) we can obtain and such that 

equation (4.4) can be reduced to 


k l ^ 1} +^ (2) ^ (3) ♦ •••) ♦ sjitf' *d 2) + # (3) ♦ •••) = - s j^ 


0 ) 


= -(s T u” 1 )r 


+j/ 3 ^ + ...) + ♦ &f' 2 ^ + * •••) = 


0) 


« -(H M K[V 


We set 


k^ 11 - -s^°> ly^.-n/ 


0) 


so that we obtain 


v n) - -K^S^ 0 * . -K L 1 s T u‘ 1 r. 


« (1) - - u ;V 0) ■ -“i'¥[ lr i 


We substitute these In equation (4.5) and obtain 

K l (s/ 2 > tJ J» ,j/ A) . ...) . S 1 (s/ 2> *fiP' * •••) 

- -s^'> ■ ♦s I u» 1 Vi 1r i 


„y 2 > .^ 3 > ^ (4 > * ...) * U v(a P> .of" >d A) * ...) 


* V 


1) 


.-U ..-1, 


.m n k l s t u v r 2 


Next put 



and obtain 


^ (2) . 
&t 2) = 


.K-^tU^nK-l 1 




r 

r 


1 

2 


Substituting Into equation (4.8) we get the revised equation 
K l (s/ 3) +jj/ 4) + ...) +■ S t (^ 3> +^ 4) + ...) = -S^ 2) 

M N (^ 3) + j / A) f ...) f U y (^ 3) +^ 4) * ...) = 
from which using the same procedure as before we can obtain 

^< 3 > . 

S< 3) . -u^h/ 2) 


(4.10) 


(4.11) 


(4.12) 


sf 

jJ 


Blrnbaum coefficients 

may be 

written as follows: 




a< 0) 

= U V r 2 


•0) 

S 


ta* 11 . 



(2) = 


a <2) . 

-UJV 1 ' 

(4.13) 

'(3) = 

-k[ 1 s 1 ^ 2) 

«< 3) . 

-u;V 2) 





It Is observed that the successive terms of the sequence can be generated 
easily by using the operators ( Sj) and (UylMfl) repeatedly on 
the result of the previous recursion. 


Truncating the Blrnbaum series at m = M*, we have (M* + 1) coeffi- 
cients for each of the rotor and stator. Considering P* points on each side 
of a blade, we have 2P* equations to determine the (M* + l) coefficients of 
the rotor and stator. Therefore, we must have 2P* = M* +1. The matrices 
Kl. Mn, Sj, and Uy are of order (2P* X 2P*). The matrices n and T2 are 

of order (2P* XI). Each of the vectors .aioj.ad) 

Is of order (2P* XI). The Blrnbaum coefficients may be written as 
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(4.15) 


- S^8) =j^ (0) - K[ 1 S t ^ 

# = ujjV 2 - M^aO « 0 {O) - u^m n ^ 

which clearly shows the coupled nature of the rotor-stator flow and the correc- 
tion terms to the respective Blrnbaum coefficients due to the aerodynamic 
Interference. The structure of the matrices K, L, M, and N Is given In the 
Appendix In truncated form along with that of ] , #2* ^3» an(1 ^4* Therefore, 
the matrices K|_, Mn, St, and Uv along with n and T 2 are given by the 
partitioned matrices 



In equation (4.16) we have K aB (R a ), L aB (R a ), M aB (R a ), N ttB (R a ), S aB (S a ), 
T aB (S a ), U aB (S a ), V aB (S a ), <ifia(Ra). ^2a(Ra). ^3a(So). and <jf 4 a (S a ) correspond- 
ing to the points R a and S a on the rotor and stator blade respectively with 
o = 1,2,3,...P*, 8 = m + 1 , m = 0,1 ,2,3, . . .M*. The matrices Kl. Mr, Sj, and 
Uy are primarily aerodynamic In nature while the matrices n and T 2 are 
purely geometric representing the blade section and the stage configuration. 

The Individual elements of the aerodynamic matrices K, L, M, N, S, T, U, V, 
and the geometric configuration matrices ^ *^3 and ^4 have been 
defined In the Appendix. 

From the simple recurrence procedure outlined above It Is possible to 
obtain the Blrnbaum coefficients of the rotor-stator combination to any desired 
approximation. In the next section we shall outline a direct method for calcu- 
lating the Blrnbaum coefficients for the rotor-stator combination directly. 


5. THE DIRECT METHOD 


In this section we shall obtain the Blrnbaum coefficients of the rotor- 
stator combination directly by satisfying the boundary conditions on the blade 
surfaces of the rotor and the stator simultaneously without resorting to recur- 
sion. For this we rewrite equation (4.4) as 


K L s/ + $ T #= 

+ Uy0 = f2 


(5.1) 


which form a pair of nonhomogeneous simultaneous equations for sf and 0. 
Solving the two equations simultaneously, jrf and 0 are given by 
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(5.2) 


* ■ ( s x\ - s\) ( s Ai - u v' r 2 ) 

a - (k['s t - M-’u v )'' (K-'r, - m n V 2 ) 

The Blrnbaura coefficient vectors si and £jjf given by equation (5.2) per- 
tain to the rotor-stator combination. Equation (5.2) gives directly an exact 
form for the Blrnbaum coefficients si and 68 given approximately In 
equation (4.15). We can relate the coefficients obtained by the two methods. 
For this we rewrite equation (5.2) using equation (3.5) as 

■ ( s ;\ -• s? *)’ 1 ( s i’ r i 

a -(V « H \) - " N V 2 ) 

where S£ and Si* are defined by 

se* = u; ] M N 


(5.3) 


(5.4) 


We can write equation (4.15) In a form similar to equation (5.3) by decou- 
pling the two equations. Thus, replacing 68 In equation (4.15a) by the second 
equation and solving for si and similarly replacing si In equation (4.15b) 
using the first of equation (4.15) and solving for 68 we get 



Comparing equations (5.3) and (5.5) we have the following relation 


( 5 . 5 ) 


(I -ifif*)"V* C ^ 1 -i ?*)' 1 


(5.6) 


between the operators of the recursive method and the direct method. Thus, 
equation (5.5) may be used to obtain the exact Blrnbaum coefficients si and 
68 of the rotor-stator combination using the basic coefficients si(Q) and 
£8(0) of the Isolated rotor and stator. Using the results obtained above. It 
Is possible to estimate the mutual aerodynamic Interference between the rotor 
and the stator. This will be described In the next section. 


6. ESTIMATION OF MUTUAL INTERFERENCE 

In the above, we have calculated two sets of Blrnbaum coefficients j/(0), 
and si, 68- As mentioned earlier, the coefficients *j/(°) , £2f(°) pertain 
to the Isolated rotor and stator and have no aerodynamic Interference effects. 
On the other hand, the coefficients si , 68 pertain to the combined rotor-stator 
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configuration In which aerodynamic Interference effects are present. All the 
flow field parameters of Interest can be calculated using the Blrnbaum coeffi- 
cients appropriate to the configuration considered. We shall define the 
Interference In terms of the axial velocity differences between the two config- 
urations. The perturbation velocity at any point can be written as the sum of 
the perturbation velocities of the rotor and stator. Assuming for simplicity 
o>P = 0 and using equation (4.1), the perturbation velocity at any point 
(r-|, e, zi) of the combined flow field can be written 


f 


^(r 1 ,e,z 1 ,t 1 ) = V x + ^ = 


A 


X 0 cos e. - X- sin e h 
2 h.j 3 n 1 


v x 


la/+ e 


loo t 
S 


x 2 sin 8|)i * X 3 cos e 
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\ 


Y 0 cos e. - Y 0 sin 0 . 

2 h 2 3 h^ 


'1 


\V 2 sin * v 3 cos 8 h; 


J 
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x» cos e. - X 0 sin e. 

2 h 1 3 h.j 


v x 
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X„ sin e. + X. cos e. 

2 h.j 3 h^ 


( 


ii 0) * 


■J 

\ 


\ 


Y- cos e. - Y_ sin 
2 h 2 3 h 2 


Y„ sin e. + Y cos e. . 

2 h 2 3 h 2 y 


V ¥ 


U 


v lot.. 
°>e s 1 


X 0 cos e. - X„ cos e. 

2 h 1 3 h^ 


V 


X. sin 0. X„ sin 0. , 

2 h^ 3 h 1 / 


( 


yta- e 


1 “s t 1 


\ 


Y_ cos 0 . - Y_ sin 0. 

2 h 2 3 n 2 


.Y sin 0. + Y cos 0 , 

V h 2 3 V 




( 6 . 1 ) 


Denoting the four terms of equation (6.1) by ^Vo* f'so* ^1s» ^1r* the P er_ 
turbatlon velocity vector V\ can be written as 


^(ri ,0,zi ,ti) = ^-o + *so ♦ ^sl + ^1r 


( 6 . 2 ) 


In equat1on(6.2) , the term ^ ro Is the perturbation velocity due to the 
Isolated rotor; Is the perturbation velocity due to the Isolated stator; 

is the Interference velocity produced by the stator on the Isolated rotor 
field; and V\x Is the Interference velocity produced by the rotor on the Iso- 
lated stator field. The total Interference velocity V\ may be written 

f\ = ^1s + ^1r (6.3) 


The Interference factor e may be defined as the ratio of the axial com- 
ponent of the combined field to the sum of the axial components of the isolated 
rotor and stator at the same point In the flow. Thus, we write e as 
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c = 


(6.4) 


^foz + ^soz 




*foz + « 


- 1 


soz 


As expected, from this equation we observe that the Interference factor 
e = (r-j, e, z -\ , t-j) Is a point function which Is periodic In time. 


7. DISCUSSION 

In the recurrence scheme given above, It Is possible to Identify the aero- 
dynamic Interference effects clearly as opposed to the direct method. Never- 
theless, the direct method Is also an efficient method for determining the 
Blrnbaum coefficients and thereby the aerodynamic characteristics. Both the 
methods discussed above enable us to calculate the overall aerodynamic perform- 
ance of the stage for a given axial and rotational velocity and a stage config- 
uration. The calculation process may be repeated for changes In the velocity 
and configuration parameters. Thus, It Is possible to calculate the perturba- 
tion changes for different stator blade angular positions <p s for a given 
rotor blade angle <p r . Similar calculations for other configuration changes 
between the rotor and stator may also be made. From a knowledge of the 
Blrnbaum coefficients, It Is possible to calculate the distribution of pertur- 
bation pressure AP r and AP S on the upper and lower surface of the blades 
containing the the aerodynamic Interference effects 


A P, 


- Yin*' * 


H cot \ > 


6 P c = 


yH 2 (^ 


% 


cot ) 


(7.1) 


We define the dimensionless axial and tangential force coefficients Cp ^ 
and Of© and the power coefficient Cp as 

c Fz ’ V<2 '- w a c > 

C FS ■ < 7 - 2 > 

C P ■ p/ <2 0*a R2 > 

The axial force coefficients CfR* and Cp$z P er un1t blade span per 
blade of the rotor and stator are given by 
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(7.3) 


dC 


FRz 


dr 


TM 


AP U Sln(a r - c ru ) - AP 4 Sln(a r - cA sin 


udu 


Jo 


*-*(>■ 


s1n(a $ - e su ) - AP^ Sln(a $ - e^)) Sin wdu> 


while the tangential force coefficients Cpp© and Cp$© per unit blade length 
per blade are given by 


dC 


dC 


f ■ - $ l h ■ 

T-i[h 


,) - Ap # COS(a - c_ # )( sin wdu 


(7.4) 


c su ) - Ap l cos(a s - "si* 


sin u dcj 


The power coefficient Cp of the rotor Is given by 
-1 /•» 


- - — f r 

R *^ 2 Jhr Jo 


AP 


u cos(a r - e ru ) - Ap^ cos(a r - c rjL ) I sin Z du dr 1 


(7.5) 


c ru = tan ~ T ul 


c rt “ tan X L1 


e su = tan_1 t U2 


£ sl = ^ an_ ^ t L2 


The upwash velocity of the rotor-stator combination Is given by the velocity 
component and may be written as 


n * ni * % 2 - * 3 ♦ v 3 a 


(7.6) 


The circulation distribution over the blade span can be obtained by Integrating 
the pressure normal to the blade chord and Is given by 


r o = - o I ( AP U - *Vr s1n “ d “ 


yM 


I 

f 


(7.7) 


F S = - ^2 | (AP U - Af Vs Sln “ d “ 
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8. CONCLUSIONS AND SUMMARY 


Using an acceleration potential method and the lifting surface theory, two 
methods, a recursive and a direct method, for solving the mutually Interacting 
aerodynamic problem of rotor-stator combination have been outlined for deter- 
mining the blade loading and other characteristics of an axial compressor stage 
with the rotor placed ahead of the stator and each with a finite number of 
blades. Expressions have been given for the resultant pressure field of the 
rotor and stator separately for an unknown distribution of surface pressure 
loading on the blades. The blade surface loading has been represented by a 
modified Blrnbaum-Glauert series with unknown coefficients. In the recursive 
method, a set of basic Blrnbaum coefficients Is obtained by satisfying the con- 
dition of flow tangency on both sides of the rotor and stator blades Individu- 
ally. The successive values of these coefficients are obtained to any desired 
accuracy by a systematic matrix procedure on the basic set of coefficients to 
account for the effects of rotor-stator Interference. In the direct method, 
on the other hand, the final values of the Blrnbaum coefficients are determined 
by a single matrix operation satisfying all the surface flow tangency condi- 
tions simultaneously. The relation between the two methods has been shown by 
relating the corresponding matrix operators. The Blrnbaum coefficients s/ 
and QH of the rotor-stator system have been related through the matrix opera- 
tors to the basic Blrnbaum coefficients j/(0) and of the Isolated 

rotor and stator. Expressions have been given for calculating the perturbation 
velocities, the circulation distribution along the blades, the axial and tan- 
gential forces on the blades and the power required by the rotor , using these 
coefficients. 
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NOMENCLATURE 


A 

a« 

B 

C 

C 

0 

hp , hj 

J k (r) 

kp , kj 

L< b > ( r ) 
a 

M 

M 

M 2 

M* 

Pr.Ps 

£ Ir( rl » ol xl) 
ds( r2» o2 t2) 

R* = M/M 

R sr * r ts/ r tr 
r(r,e,z) 

H * r/r tr 
r 2 - v Tkt r l 
r hr* r hs 

r tr* r ts 

t 


arbitrary constant of Integration equation (3.1) 
free stream speed of sound 

arbitrary constant of Integration equation (3.1) 
blade chord 

arbitrary constant of Integration equation (3.1) 

arbitrary constant of Integration equation (3.1) 

rotor, and stator (hub/tip) radius ratio 

Bessel function of the first kind equation (2.3) 

rotor and stator circumferential mode numbers 

generalized Laguerre function, of degree a and order b 
equation (2.3) 

Mach number of axial flow velocity = W a /a,„ 

rotor tip Mach number = Cr^p/a,*, 

Mach number of swirl component at rotor exit 

number of Blrnbaum series coefficients considered 

radial pressure functions of rotor and stator 

perturbation velocities of the rotor In the helical coordinate 
system 

perturbation velocities of the stator In the helical coordinate 
system 

characteristic radius of rotor 

stator tip radlus/rotor tip radius 

position vector of a point In cylindrical coordinates 

dimensionless radial coordinate 

radius parameter defined In equation (2.3) 

hub radius of rotor and stator 

tip radius of rotor and stator 

time 
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t*| = t/to 

*o * r tr/'" l a 
^u,v,w) 

k&rKK) 

V r ,V s 

(^Vr.^Ve^rz) , 
(^sr.^se.^sz) 

W r ,W s 

Yk(r) 

z ro» z so 

z rl > z r2 
z sl * z r2 
z r* z s 

Z 1 = z / r tr 

z' 7* 

4 ur »z us 

z Lr* z Ls 

z cr * z cs 

z Tr« z Ts 

A A 

^r.^s 
a r » a s 
a 2r 

B = (1 - M2)l/2 

b 2 = (1 - nW /2 

APr »Ap s 

c 

Y 


dimensionless time 
characteristic time parameter 

local velocity vector In cylindrical coordinates (r,e,z) 
perturbation velocities of rotor 
perturbation velocities of stator 
circumferential velocity of fluid for rotor and stator 

(r,e,z) Components of perturbation velocity of rotor and 
stator 

axial velocity of fluid for rotor and stator 
Neumann function of order k 

axial position of mid-rotor/stator plane from reference 
origin 

axial coordinate of rotor leading and trailing edge 
axial coordinate of stator leading and trailing edge 
number of blades In rotor and stator 
dimensionless z coordinate 

upper surface ordinate of rotor and stator blade profile 

lower surface ordinate of rotor and stator blade profile 

mean line ordinate of rotor and stator blade profile 

local half thickness of rotor and stator blade profile 

parameter defined In equation (2.6) 

stagger angle of rotor and stator blades 

exit blade angle of rotor 

Prandtl parameter 

parameter used In equation (2.5) 

net pressure difference on rotor and stator 

rotor-stator Interference factor 

ratio of specific heats of gas 
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e 

r r* r s 

<p r ,(p$ 

<Pr ><Ps 
*kl( r ) 

$* 

nkl(r) 

Pr*P$ 
w r t o) s 

Hi 

A kt 

nki 

£r>fs 

rsj /v 

u r ,w s 


G 

Subscripts: 

Re/lm 

r,s 

h,t 

l.j.k 

r,e,z 

r,a,T 


Azimuth angle of field point In cylindrical coordinate system 

strength of local rotor/stator bound vortex 

azimuth angle of pressure pole/dipole on rotor, stator blade 

mean offset angle of first rotor and stator blade 

normalized rotor eigenfunction of mode number k and P-th 
eigenvalue 

rotor subsonic pressure function equation (4.10) 

normalized stator eigenfunction of order k, degree T 
and l-th eigenvalue 

stator subsonic pressure function 

radial position of pressure pole/dipole on rotor, stator 
blade 

Glauert angle of blade defined In equation (2.9) 

l-th radial eigenvalue of rotor of circumferential mode 
number k 

l-th rotor axial eigenvalue of mode number k 
stator radial eigenvalue 

axial coordinate of rotor/stator pressure pole 

angular frequency of oscillating rotor/stator pressure pole 
(rad/sec) 

angular velocity of rotor 

real/imaginary part 
rotor, stator 
hub, tip 

unit vectors along (x,y,z) directions 
cylindrical coordinates 
local helical coordinates 
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APPENDIX 


The functions ^ ^ In equation ( 3 . 2 ) and the coefficients 

K r and K s contained In It are defined by the following equation (A-l) 


f MH , kl c „ n : 

'o-^^T s « n Z r 


, ■£' .. 


1T 1 

f 2 e l dZ 1 


f 0 = — - f rt tan 
2 p r o 2 r 

k r f Z l IT, 
f 3 = R*J f 2 e dZ l 

•'m 


- z i/ s 


°l = z, /sin 6 h , 


(t 


*2 MrZ C0S e h2 * s1n 9 h2 Z 1 


2 r ■ 2 ! - <r 


T o ■ - T tan “ 


r exp 1 |» r C r * k r (e ' - * r - 57 | 


T 1 * Z l( f o * “r * «;)- Vr * k r (9 ro ' *r ' *r> 

( 9 so * " is "' s ) ' “ s N * “ s ° 2 S ' n < ’ 1 ' 2 


T 2 ' k s 9 so * 
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T 3 = k s l 6 ' *s ‘ *s ‘ R** ) + "s Sl ° 2 0 h2 l 1 * R^* C0t 0 h2j f s 


f 1 * 5^7 cot vj* 


tan e hl = H/flr^ 
= ^O g l 


% = *^ 40 9 4 


d ^0 X kl 


\A*[ J. , (X. 0 r, ) - 


dr i ■ Vki '•* k - ] kn 


* , k+l (x kl r l )1 + Y k-l (x ka r l ) ‘ Y k+l (x kl r l ) 


•?0 = *ki< r l>Ai 

- -1/2 r, 

, «5o ± : 

4° - dr, - n ?kp 


(h^ - l) L i 2i) < r 2> - Hf L lO >(r 2>] 


+ r. 
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(V’W' 
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f 2 e l dZ 1 


1 cos 26 
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6 ro = 6 - R* 


Z s -" z i - <s 


tan e h2 = M ^ r l ~ M 2^ = ( cot e hl ~ tan a 2r^ 


-1 


*2 = * 20 q 2 
*5 = *50 9 5 

*50 s *lkt^ r l ) r< *2 e 


f 5 ■ T % * ^ tan a s S9n h 


■, • r 

J oo 


IT, 


f c sgn Z e dZ, 
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IV 2 e ' '“r 9 l) 
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' 00 
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1T 2 

f 5 e z d Zl 
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The matrices K, L, M, N, S, T, U, and V are expressed by a linear 
combination of the functions X-j , X 2 , X 3 , Y-| , Y 2 , and Y 3 using 
the coefficients A^j. The matrices K, L, M, N, S, T, U, V, and 
coefficients A^j are given In equations (A-2) and (A-3) as follows: 

K = f 1 = A n X 1 + fi 12 X 2 ♦ *-, 3 X 3 S - g 1 = S 14 Y 1 + A 15 Y 2 + S 16 Y & 


*21 *1 + *22*2 * *23*3 


*31*1 + *32*2 * *33*3 


*24 Y 1 * *25 Y 2 + *26 Y 6 


*34 Y 1 * *35 Y 2 * *36 Y 6 


f 4 “ *41*1 * *42*2 + *44*4 


9 4 * * 44 Y l + *45 Y 2 + *46 Y 6 


A^ = COS 0 ( t yi sin a f + COS a f ) 

& 12 = sin e hl ( T U1 cos a r " s * n a r^ " s1n 6 cos e hl^ T ul s1n a r + cos a r^ 

5 13 = cos 9^1 ( T ui cos a r " sln a r ) + s1n ® s1n e hi^ T ui s ^ n a r * cos a r^ 

5 14 = cos e (t^ sin <* r ♦ cos a f ) 

*15 = s1n °h2 (x Ul cos a r “ s1n a r^ ~ s1n e cos e h2^ T Ul s1n °r * cos “r* 

A 16 = cos ©h2^ T Ul cos a r “ s ^ n °r^ * s1n ® s ^ n e h2^ T Ul s1n a r + cos a r^ 

*<f, = sin 0 cot 0 k1 (T. n sin <* + cos a ) - (t m1 cos a - sin a ) 

1 hi ui r r ui r r 

A 3 i = cos ®( t u 2 “$ + cos 

A 32 = sin ®^,i ( T u2 cos a s “ sln a s^ ' s ^ n e cos e hl^ T U2 s * n a s + cos a s^ 

A 33 = cos ®f,i ( T (j 2 cos a s " s ' n a s^ + s1n ® sln e hi^ T U2 S ^ n a s + C0S a s^ 

A 34 = COS 0(t U 2 sin a $ + COS a $ ) 

A 35 = sin ®^2 ^ t U2 cos a s " sln a s^ ” s1n e cos e h2^ T U2 s1n a s * cos a s^ 

*36 = cos e h2^ T U2 C0S a s " s ' n °s> + s1n 6 s1n 6 h2^ T U2 s1n a s f C0S a s^ 

^ = sin 6 cot %2 ^ t U 2 s ^ n tt s * C0S a s^ " ^ T U2 cos a s " s1n a s^ 
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cos e (t l1 sin a f + cos a r ) 

sin e hl (T U cos a p - sin cy) - sin e cos e hl (T u sin cy + cos oy) 
cos e hl (T u cos a r - sin cy) + sin e sin ® hl (-r L1 sin <y + cos a r ) 
COS e (t l1 sin a r + COS a r ) 

sin 6 h2 ^ T Ll cos a r " s1n a r^ " sln 6 C0S e h2^ T Ll sin a r + cos a r^ 
cos © h2 ( T L] cos a r " s ^ n a r ) + s1n 6 s1n e h2^ T Ll s1n a r * cos a r^ 
sin e cot © hl ( T L1 s1n a r ♦ cos a r ) ~ ( t li cos a r ” s1n a r^ 
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sin cy ■ 
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COS a 
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6(t L2 
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V T 
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V T 
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^ - ill! ) - 5 111 SF IU5 °hP T L2 $ 


hi '12 


^ - din at^ ) t din v din L2 S i 


sin e cot e k 2 ( T L2 s1n a s + cos a s ) - ( t l 2 cos a s “ s ^ n a P 


The schematic structure of the matrices K, L, ...U, and V for 
determining the (M* +1) Blrnbaum coefficients of the rotor and stator Is 
shown In equation (A-4) below: 
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